Tight-binding parametrization for the chromium nitride: A NMTO study 
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We investigate the band structure of chromium nitride using the iVth-order muffin-tin orbital 
(NMTO) based downfolding technique. The effective hopping Hamiltonian parameters are obtained 
using NMTO downfolded basis sets, which consist of Cr d bands including e 9 and t% g states for both 
cubic and orthorhombic lattice. We analyze the chemical bonding and tight-binding parameters 
from the tight-binding Hamiltonian, further the effect of lattice distortion is discussed according to 
these parameters. 
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I. INTRODUCTION 



Chromium nitride (CrN) has great potential for indus- 
trial applications as hard, wear-, and corrosion-resistant 
coatings. It also attractes much attention due to its in- 
teresting fundamental physical properties and nonuni- 
form picture of the electronic structure 1 - - — . The mag- 
netic and structural properties of CrN have been well de- 
fined: at room temperature CrN is paramagnetic (PM) 
in the rock-salt (RS) structure (Space group Fm3m), 
but below Neel temperature T/v~286K it becomes an 
antiferromagnetism(AFM) orthorhombic phase (Space 
group Pnma)^-. The AFM magnetic ordering consists 
of ferromagnetic(FM) planes, stacked antiferromagnet- 
ically along the [110] direction every two FM Cr lay- 
ers (AFM? 110 ,), which has been identified by neutron 

scattering^. In contrast to the general agreement on 
the magnetic ordering and crystal structure across Tjv, 
the reports on its electronic properties are quite con- 
tradictory, with (i) metal to metal^i 7 -, (ii) insulator to 
insulator 3 *!!, and also (hi) insulator to metal transition 8 . 
A recent study has identified a softening of bulk CrN 
under pressure due to a manifestation of a strong com- 
petition between different types of chemical bond that 
are found at a crossover from a localized to a molecular- 
orbital electronic transition!. Moreover, experiment also 
indicates that the N-vacancy concentration and crys- 
talline defects strongly affect electron transport^. There- 
fore the stoichiometric CrN has been considered a cubic 
paramagnetic correlated insulator at room temperature, 
but an orthorhombic antifcrromagnet metal below Tn&z.. 

However, on the theoretical side, First-principles den- 
sity function theory (DFT) calculation a 12 ' 13 predict 
metallic ground states for both cubic and orthorhombic 
lattice. The results show Cr 3d partial density of states 
(DOS) at the Fermi level (Ep), but no gap observed, 
which is in contrast to the experimental measurement 4 . 
After considering the On-site Coulomb repulsion for Cr 
3d electrons, DFT+U calculations for the AFM phase 
could open a direct gap when U value is larger than 3 
eV— . Moreover, in agreement with experimental obser- 
vation, the DFT calculations 1 ^ have confirmed that the 



ground state of orthorhombic CrN is AFM? 11Q , . It also 
predicts the ground state of the cubic lattice is antifer- 
romagetism in which every one FM Cr layer stackes an- 
tiferromagnetically along the [110] direction (AFMf 110 ,). 
However, in contrast to the theoretical results, the exper- 
iments show that cubic CrN remains a paramagnetic in- 
sulator over the entire measured temperature range of 10- 
295 K 3 . Therefore, the relationship between magnetism 
and structure as well as the related electronic structure 
of CrN is still a controversial issue. In this paper, we em- 
ploy the recently developed ./Vth-order muffin-tin orbital 
(NMTO) method 1 ^— to produce real space Hamiltonian 
parameters for both cubic and orthorhombic phases, fur- 
thermore investigate the hopping integrals of the tight- 
binding Cr d bands and analyze the chemical bonding. 



II. THE CRYSTAL STRUCTURES AND 
COMPUTATIONAL METHOD 

The two kinds of structures are very well known and 
are described as bellow: One is the cubic rock-salt struc- 
ture with a lattice constant of a — 4.148 A 1 . Each Cr 
atom is surrounded by six N atoms, thereby providing the 
octahedral environment at the Cr site, which leads to the 
splitting of the degenerate d orbital into t 2g and e g states. 
Every Cr atom has 12 first nearest neighbor (first-NN) 
Cr atoms and 6 second nearest neighbor (second-NN) 
Cr atoms. Below the Neel temperature, CrN has the 
orthorhombic antiferromagnetic structure (lattice con- 
stants are a = 5.757 A, b = 2.964 A, c = 4.134_A)!£, 
in which the lattice expands about 1.8% along [110] di- 
rection and shrinks about -1.1% along [110] direction of 
the original cubic structure, respectively. Since there is 
the simultaneous occurrence of both magnetic order and 
structure transition, in order to separate the effect of 
magnetism and structure distortion, we ignore the mag- 
netic ordering for the orthorhombic phase. In the or- 
thorhombic lattice, each Cr atom is still sixfold coordi- 
nated by nitrogen atoms and the corresponding octahe- 
dra is distorted due to the tiny shear deformation in xy 
plane, while c-axis is roughly unchanged. The 12 first- 
NN Cr-Cr pairs are split into two longer and two shorter 



pairs in xy plane and roughly unchanged 8 Cr-Cr pairs 
out of xy plane. It is also noted that the bond distance 
of Cr-N as well as second-NN Cr-Cr pairs remain nearly 
constant. 

In the present study we firstly use the tight-binding 
linear muffin-tin orbital (TB-LMTO) method with the 
atomic-sphere approximation techniqu e) 15 ' 16 to generate 
electronic structure. The non-spin-polarized calcula- 
tions for both the cubic and orthorhombic CrN phases 
have been performed within Local Density Approxima- 
tion (LDA). Then NMTO method is used to carry out 
the localized Wannier function which is the orthogonal- 
ized Lowdin functions from the atomic orbitals. The 
downfolding technique in the NMTO method allows to 
produce minimal bands which follows exactly the bands 
derived with the large basis set, which actually is the 
maximum localized Wannier function. We choose all of 
Cr d orbitals to form the minimal basis set, and further 
the corresponding hopping integrals of the tight-binding 
Cr d bands arc derived from the NMTO calculations. 



III. RESULTS AND DISCUSSIONS 
A. LDA band structure 

The non-spin-polarized band structures of CrN are cal- 
culated with the full atomic orbitals, which are shown in 
the left-panels of Fig.l and Fig. 2 for the cubic lattice 
and the orthorhombic lattice, respectively. In the cu- 
bic lattice, the high-symmetry points in Billouin Zone 
are: L = (0.5,0.5,0.5), G = (0,0,0), X = (0,1,0), 
W = (0.5,1,0), L = (0.5, 0.5,0.5),K = (0,0.75,0.75), 
G = (0,0,0) and Z = (0,0,1). In the orthorhombic 
lattice, the high symmetry points in Billouin Zone are: 
L = (0.5,0,0.5), G = (0,0,0), X = (0.5,0.5,0), W = 
(0.75,0.25,0), L = (0.5,0,0.5), K = (0.375,0.375,0.75), 
G = (0,0,0), and Z = (0,0,1). The chosen high sym- 
metry lines would pass through similar way in Billouin 
Zone of both the cubic and orthorhombic lattices, so that 
we could compare the discrepancy of band structures be- 
tween two crystal structures. 

The electronic structure of CrN is extensively inves- 
tigated previously. The strong hybridization has been 
found between Cr d and N p states. The wide band 
indicates the strong p-d bonding between Cr d orbitals 
and N p orbitals. Comparing the bands shown in Fig.l 
and Fig. 2, it is found that the band structures of both 
phases are quite similar except that the small band split- 
ting about 10 meV induced by orthorhombic distortion. 
As shown in Fig.l, the lowest three bands are primarily 
p states of N distributing from -9 eV to -2 eV. The t 2g - 
derived states at the Cr site span the energy of about -3.9 
eV^0.9 cV, while the e s -derived states of Cr occupy the 
energy range from -0.7 eV to 4.2 eV. Around the Fermi 
level, the bands are primarily Cr 3d states of t 2g mani- 
fold, with e g bands situating in higher energy due to the 
crystal-field splitting. In the cubic lattice, the d-orbitals 



are simply split to t 2g (d xy , d yz , d xz ) and e g (d 3z 2_ 1 , 
d x 2_ y 2) mainfolds due to symmetry of the CrNg octahe- 
dra. While in the orthorhombic lattice, as the shear de- 
formation in xy plane reduces the symmetry of distorted 
octahedra, t2 g is further split to d xy , d yz and d xz states, 
and e g is split to d^-i and d x 2_ y 2 states, respectively. 



B. Downfolding onto the Cr-d manifolds for the 
cubic lattice 



Since the band gap observed experimentally in para- 
magnetic cubic phase could not be captured by LDA or 
GGA calculations, it suggestes that CrN should be elec- 
tron correlated system^. Therefore, constructing a model 
hamiltonian which could be used in many body picture 
to perform the realistic calculation is highly expected. 
We propose a downfolding scheme in which all orbitals 
of whole atoms except Cr-d are downfolded. Both Cr t 2g 
and e g states are considered to constitute the effective 
orbitals because electrons occupation on higher e g state 
could not be ignored in both cubic and distorted struc- 
tures. For the local coordination, the x and y coordinate 
axis of CrNg octahedra are chosen in original direction 
of the coordinate axis of the cubic lattice. The effective 
hopping Hamiltonian matrix of dimension 5x5, built 
up by the five d effective orbitals. Fourier transforma- 
tion of the downfolded Hamiltonian (H(k)— >H(R)) gives 
a TB electronic Hamiltonian in real space consisting of 
hopping over up to five NNs. For many studies, it is de- 
sirable to have a shorter-range Hamiltonian which can be 
achieved by the downfolding technique. Taking the first 
nearest neighbors and the second nearest neighbors into 
account, the total Hamiltonian matrix can be expressed 
as follow: 



H = £ 
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In this formula, R n and R n are the coordinates of 12 
first nearest neighbors and 6 second neighbors, respec- 
tively. Here only one representative hopping matrix at 
each NN is shown as bellow, but due to the crystal sym- 
metry all the other NN hopping integrals can be derived 
from proper unitary transformation. 
The basis set of Cr-d NMTO orbitals: 

|X X ) = {\xy), \yz), \zx), |3z 2 - 1), \x 2 - y 2 )} (2) 

The on-site hopping integral (The unit is eV and Ep 
= 3.2685eV.): 



/ 3.1618 
3.1618 



,000 



















3.1618 
5.3497 











\ 







5.3497 / 



(3) 



The first nearest neighbor (The unit is meV): 
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The second nearest neighbor (The unit is meV): 
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The on-site term t^, m is diagonal, and doubly de- 
generated e g and triply ti g are presented. The crystal- 
hied splitting between e g and t2 g state is 2.2 eV. In the 
downfolding procedure, we have obtained the hopping 
matrix elements of all Cr-d orbitals, thus we can derive 
the energy bands from these hopping matrixes through- 
out the Brillouin zone. Hopping integrals further than 
third NN are small. The largest values are 35 meV for 
third NN, 59 meV for fourth NN and 54 meV for fifth 
NN, respectively. The hopping of third NN should be 
pass through chains like Cr[000]-Cr[||0]-Cr[110] and 
the hopping of forth-NN should be pass through chains 
like Cr[000]---Cr[±±0]---Cr[l±±], so that they are small. 
The others of matrix elements are mostly zero. Further 
NN hopping integrals are negligible. 

From the effective TB second-NN Hamiltonian, the 
best optimized downfolding bands are plotted with black 
dash lines in the right panel of Fig.l. The bands obtained 
by downfolding all the other channels expect the Cr-d 
channels include only the first- NN and second-NN Hamil- 
tonian in real space. There is distinguish discrepancy 
between the full bands and downfolding bands, which in- 
dicates that the long-range hopping have certain contri- 
bution. We find that the third NN and fourth NN terms 
are necessary to fit the full bands. 

Furthermore, the second-NN Hamiltonian derived 
from downfolded bands are expressed in terms of the 
Slater-Koster integrals dda, ddir, ddS, which could be 
used to understand the chemical bonding between Cr 
atoms on different sites. Starting from the first-NN 
hopping parameters of t2 g , wc obtain ddcr=532meV, 
rfd7r=4meV, and dd<5=180meV, while from the first-NN 
hopping parameters of e g we obtain dd7r=168meV and 
jdda + jddd =21meV . The difference of chemical bonding 
obtained between ti g and e g orbitals is due to the influ- 
ence of N-2p obitals. Because the e„ orbitals have strong 
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FIG. 1: (Color online) The band structure calculated with full 
basis (left panel) and the downfolded Cr-d bands (right panel) 
of cubic CrN. In the right panel also shows the best-fit TB 
bands with hopping integrals up to second-NN (black dashed 
lines). The LDA d bands are also shown for compasion (red 
solid lines). 



ey<7=-530meV, <i<i<5=32.7meV for e g orbital, respectively. 
The chemical bonding between e g orbitals are dominated 
by dda plus a little dd5 without ddir, while between ti g 
orbitals are ddir dominated. 

It could be found that the effective bonding connected 
to second-NN Cr atoms are actually indirect bonding, 
while for first-NN Cr atoms are mainly direct bonding 
between Cr d orbitals, if we perform downfolding of the 
Hamiltonian in the larger basis including N 2p orbitals. 
The second-NN indirect bonding between ti g orbitals are 
primarily dominated by the nitrogen-mediated pdn cou- 
pling, while bonding between e g orbitals are primarily 
dominated by the nitrogen-mediated pda coupling, along 
x, y and z-axis directions via the Cr-O-Cr chain. As 
the pda and pdir bonding are very strong (~ 1.5 eV), 



so that the indirect bonding (~ * d ) of second-NN has 
comparable strength with the direct dd bonding of first- 
NN. 



Downfolding onto the Cr-d manifolds for the 
orthorhombic lattice 



We now turn to the orthorhombic lattice. As same as 
the cubic case, the effective TB Hamiltonian is of dimen- 
sion 5x5, which is defined on the basis of the effective 



pda antibonding with N-2p states, while t<2 g orbitals form t<2 g and e g orbitals for the Cr sites of the orthorhombic 



weak antibonding pdn coupling with N-2p states. Since 
the influence of N-2p tail is captured in the NMTO down- 
folding procedure, therefore the two-center integral treat- 
ment of Slater-Koster integrals formation provide differ- 
ent bonding strength. Moreover, from the second-NN 
hopping parameters, we obtain that the chemical bond- 
ing are dd7r=-117meV, ddS—9.8meV for tig orbital and 



lattice. The hopping integrals of Cr-d Wannier orbitals 
in the real space are presented up to the second NN. 
Since the Cr atoms centered CrN 6 octahedra is slightly 
distorted by the orthorhombic distortion, the 12 first-NN 
terms in Hamiltonian are split into three groups accord- 
ing to crystal symmetry: two hopping integrals in xy 
plane connect to (55O) and (-5-5O) sites, other two hop- 
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The second nearest neighbor (The unit is meV): 



FIG. 2: (Color online) The band structures calculated with 
full basis (left panel) and the downfolded Cr-d bands (right 
panel) of orthorhombic CrN (Pnma). The right panel show 
the best-fit TB bands with hopping integrals up to second- 
NN (black dashed lines). The LDA d bands are also shown 
for compasion (red solid lines). 



ping integrals in xy plane connect to (5-5O) and (-55O) 
sites, and the rest eight hopping integrals are correspond- 
ing to the sites out of xy plane. Similar with cubic case, 
only one representative hopping matrix for each group is 
shown in the following. 

For convenient, the a and b axis are further rotated 45° 
in NMTO calculations, pointed to (110) and (110) of the 
orthorhombic lattice, while the c axis is unchanged. In 
this local coordination system, the a, b and c axis are as 
same as that in the cubic structure, thus we can compare 
the hopping integral directly. 

The basis set of Cr-d NMTO orbitals: 

\ X ^)^{\xy). / \yz),\zx) : \3z 2 -l) : \x 2 -y 2 )} (6) 

The on-site term (The unit is eV and E F = 3.4286 
eV.): 
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The first nearest neighbor (The unit is meV): 
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The best optimized downfolded bands are also plotted 
in the right-hand panel of Fig. 2 (black dash lines). Sim- 
ilar with that of the cubic lattice, we could find the dis- 
crepancy between the full bands and downfolded bands, 
because the bands obtained by downfolding only take ac- 
count of the first-NN and second-NN Hamiltonian in real 
space. 

It is found that the on-site terms are not diagonal un- 
less in { 0.9998d X8/ +0.0158d 32 2_ 1 , ^d xz +^d yz , &d xz - 

^d yz , -0. 0158d xy +0. 9998d 3z 2_ 1 , d x i_ y -i } representa- 
tion, corresponding to eigenvalue {3.3106, 3.3188, 3.3266, 
5.5604, 5.5641}, respectively. Since the orthorhom- 



bic distortion in the xy plane, the d xz and d yz or- 
bitals are not orthogonal, while the orbital combinations 

2^ xz ~^2^v z an d 2 ^z- 2 dyz ar e orthogonal, whose 
directions point to the a and b axis of the orthorhom- 
bic lattice, respectively. The orbital combinations 
of 0.m§M xy +0M'oM3 Z 2_ 1 and -0.0158d a; y+0.9998d 32 2_ 1 
indicate that the pure state of d xy or d 3z 2_ 1 orbital in the 
cubic lattice slightly mix each other due to the distortion 
effect induced by the expansion of a axis and shrinking 
of b axis in the orthorhombic lattice. 

The influence of the orthorhombic distortion on the 

hopping integrals of first-NN and second-NN is clearly 

1 io 
indicated, when we compare the hopping terms t^f m , 

- — -0 0-- 

t L',m ' t m\m and Sn'.m of tne distorted lattice with 
the corresponding hopping terms of the cubic lattice. 
It indicates that most of first-NN hopping terms have 
a approximated relationship: t^ l , m (orthorhombic) + 
*m^m (orthorhombic) ~ 2 x t^f m (cubic), while second- 
NN hopping matrix tf®} m change a little. Thus, the ef- 
fect of orthorhombic distortion on electronic structure is 
almost offset in these terms. The net contributions pri- 
marily originate from the first-NN hopping terms of t xz . yz 
and t^ z 2_i xy between the Cr atoms in xy plane, which 



also indicate the nature of the orthorhombic distorted 
effect similar with that found in on-site terms. There- 
fore, the total energy should change very tiny when the 
lattice of CrN distorts from the original cubic to the or- 
thorhombic structure, as been shown in the total energy 
results of the previous nonmagnetic first-principles elec- 
tronic structure calculations 12 . 



ing to the TB parameters. It is found that the first-NN 



hopping integral t xzyz and t 3z 2_ 



l,xy 



should be respon- 



sible for the effect of orthorhombic distortion and long 
range hopping terms have little contribution. Therefore, 
the non-spin-polarized result with Cr d orbitals involving 
second-NN hopping terms might be a good starting point 
for the many-body calculations to consider the complex 
physical properties. 



IV. SUMMARY 



In conclusion, by using the downfolding technique with 
the NMTO method, we have investigated the electronic 
structure of both the cubic and orthorhombic structure 
CrN, and obtained the parameters for the TB Hamilto- 
nian with downfolded Cr d orbitals. Further the chemical 
bonding and the distortional effect are analyzed accord- 
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